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ABSTRACT 

We explore the energetics and plasma composition in FR II sources using a new 
simple method of combining shock dynamics and radiation spectrum. The hot spots 
are identified with the reverse shocked region of jets. With the one-dimensional shock 
jump conditions taking account of the finite pressure of hot ICM, we estimate the 
rest mass and energy densities of the sum of thermal and non-thermal particles in 
hot spots. Independently, based on the Synchrotron Self-Compton (SSC) model, we 
estimate the number and energy densities of non-thermal electrons using the multi- 
frequency radiation spectrum of hot spots. We impose the condition that the obtained 
rest mass, internal energy, and number densities of non-thermal electrons should be 
lower than those of the total particles determined by shock dynamics. We apply this 
method to Cygnus A. We examine three extreme cases of pure electron-positron pair 
plasma (Case I), pure electron-proton plasma with separate thermalization (Case II), 
and pure electron-proton plasma in thermal-equilibrium (Case III). By detailed SSC 
analysis for Cygnus A and 3C123, we find that the energy density of non-thermal 
electrons is about 10 times larger than that of magnetic field. We find that the Case 
III is not acceptable because predicted photon spectra do not give a good fit to the 
observed one. We find that Case II can also be ruled out since the number density 
of non-thermal electrons exceeds that of the total number density. Hence, we find 
that only pure plasma (Case I) is acceptable among the three cases. Total kinetic 
power of jet and electron acceleration efficiency are also constrained by internal energy 
densities of non-thermal and total particles. 

Key words: Radio Galaxies: general — shocks: theory — radiation mechanisms: non- 
thermal 



1 INTRODUCTION 

Energetics and plasma composition of relativistic jets in ac- 
tive galactic nuclei (AGN) is one of the unresolved issues 
for understanding the jet physics. Recently, Chandra X-ray 
Observatory has detected inverse Compton X-rays from hot 
spots in FR II sources (e.g., Harris et al. 2000; Wilson, Young 
& Shopbell 2000; Hardcastle, Birkinshaw, & Worrall 2001; 
Hardcastle et al. 2002). We can now remove the degener- 
acy between the energy density of relativistic electrons and 
magnetic field by using both the inverse Compton compo- 
nent and the synchrotron component. 

Since these photons are produced by shock-accelerated 
electrons, the existence of thermal particles which drive the 
shocks is not directly probed. In this paper, we will explore 
the physical quantities of thermal plasma by using the shock 
dynamics. It is clear that both thermal and non-thermal 
plasma contribute to the dynamical processes in hot spots. 



Many authors have studied the dynamics of cocoon using 
both analytical and numerical approaches (e.g.,Begelman & 
Cioffi 1989; Falle 1991; Kaiser & Alexsander 1997; Smith et 
al. 1985; Scheck et al. 2002). However, little has been dis- 
cussed about the energy and number density ratio of the 
shock-accelerated electrons to the total particles (thermal 
plus non-thermal) and the issue of electron acceleration ef- 
ficiency is still open up to now. If the energy density of ac- 
celerated particles is negligible compared to that of thermal 
particles, it means that there is a large amount of missing 
kinetic power of thermal particles. Furthermore, quantita- 
tive study of missing power will give some constraints on 
the problem of plasma composition. 

In order to explore these issues, in §2, we outline the 
method to constrain the energetics and plasma contents by 
comparing the rest mass and energy densities of shock accel- 
erated electrons and total particles obtained by the analysis 
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of shock dynamics. Here we include the issue of shock jump 
conditions. In §3, we investigate the energetics of shock ac- 
celerated electrons, magnetic field and radiation field in hot 
spot of Cygnus A and 3C123 based on observed non-thermal 
radiation spectra. We apply the method developed in §2 to 
Cygnus A. We summarize our conclusions of constraints on 
energetics and plasma composition of FR II sources in §4. 



2 GENERAL CONSIDERATION 

It is widely accepted that shock waves accelerate some frac- 
tion of the electrons and protons in the shocked region and 
then shock accelerated (relativistic) electrons produce the 
observed non-thermal emission. It is worth emphasizing that 
by definition, physical quantities of shock accelerated parti- 
cles never exceed those of total particles. Based on this key 
point, the aim of this paper is to constrain the energetics 
and plasma contents. 

To this end, here we will investigate the number and 
energy densities of total particles in the shocked regions (see 
Figure using the simple shock jump conditions with the 
aid of observed ICM physical quantities and advance speed 
of hot spots. By comparing the number and energy densities 
of relativistic electrons and total particles (we will explore 
the number and energy densities of the shock accelerated 
electrons in the emission region of hot spots in the next 
section), we give constraints on the energetics and plasma 
contents. 



2.1 Shock Dynamics Applied to the Hot Spots 

Here we discuss the dynamics of shocks at the head of jets. 
Let us consider the dynamics of an AGN jet which impinges 
into a hot Intra-Cluster Medium (ICM). FigureQshows the 
schematic view of interaction between the AGN jet and ICM. 
We know that two shocks form: the forward shock that prop- 
agates into the ICM and the reverse shock that propagates 
into the jet. A contact discontinuity separates the shocked 
jet and the shocked ICM. 



2.1.1 Basic Assumptions 

In order to grasp the essence of shock structure of AGN 
jets, we make several assumptions: (1) we use the ID shock 
jump conditions, (2) we regard the forward shock as a non- 
relativistic one since the advance speed of hot spots is es- 
timated to be in the range 0.01c to 0.1c (e.g., Liu, Pooley 
& Riley 1992; however see also Georganopoulos & Kazanas 
2003), and adopt the adiabatic index of equation of state as 
5/3, (3) we assume the reverse shock as a relativistic one 
although the jet speed on Mpc scales is still open. Some 
are suggested to be relativistic (e.g., Tavecchio et al. 2000), 
(4) we assume that the magnetic fields are passive and ig- 
nore their dynamical effects, and, (5) we treat only one- 
component plasma compositions (i.e., pure or pure ep) 
for simplicity. The validity of the assumption (1) is discussed 
at the end of this section. 



2.1.2 Shock Jump Conditions 

As shown in Fig. 1, we use the terminology of region i (i=l, 
2, 3, and 4) with the number labeling the four regions in the 
head part of AGN jet as follows: (1) the unshocked ICM, 
(2) the shocked ICM, (3) the shocked jet which is identified 
with hot spots, and (4) the unshocked jet. Fluid velocity 
and Lorentz factor in the region i measured in the ICM 
rest frame are expressed as Vi — PiC, and Pi, respectively. 
Relative velocity and relative Lorentz factor between the 
region i and j (velocity of region i measured from region j) 



are expressed as Vij = 



^ PijC= -Pj. 



and Pi, = P, 



respectively. As for the position of the forward shock front 
(FS), the contact discontinuity (CD), and reverse one (RS), 
we use the same labeling (i=FS, CD, and RS). Each region 
is characterized by three physical quantities; the rest mass 
density pi, the pressure Pi, and the velocity Vi. In order to 
make the argument independent of the plasma composition 
of jets, we use the rest mass density rather than the number 
density. In the next subsection, we will discuss the issue of 
number density by specifying the composition as pure or 
ep plasma. 

Within the framework of ID planar shock, pressure and 
velocity are uniform in each shocked region. Then, along the 
CD, we have velocity and pressure balance as, V2 — «3 and 
P2 = P3. In general, we can solve for 3-1-3 = 6 physical 
quantities p2, P2 = P3, V2 = "3, Ps, wfs, and vks, when 
3 + 3 = 6 upstream quantities such as pi, Pi, ui, p4, P4, 
and Vi are given. In the case of the shock conditions of ac- 
tual FR II sources, it is convenient to choose 6 givens in a 
different way. The properties of ICM (upstream) are known 
from X-ray observations to give Pi = Picm, pi ~ Picm, and 
wi = 0. The hot spot (downstream) advance speed vus is 
inferred from observations. We regard that the unshocked 
jet is cold and P4 = 0. Although there are some amount 
of relativistic electrons in the jet, observations and numer- 
ical simulations suggest that the internal Mach number of 
Cygnus A jet should be high (Carilli & Barthel 1996). Fur- 
thermore, to obtain a conspicious hot spot identified with a 
strong shock, dynamically cold jet matter is the most nat- 
ural choice. We adopt P4 as a free parameter. To sum up, 
for actual FR II source, we regard that (1) V2 — V3 — vus 
is given (observable) and that (2) p4 can be solved. Then, 
by using the plausible shock conditions, we can obtain p2, 
P2 = P3, vfs, P3, v-Eis, and p4, as functions of P4. 

Using the well-known shock jump conditions in non- 
relativistic limit (Landau & Lifshitz 1959), the forward 
shocked region quantities are given as 



vfs 



P2 



P2 = 



+ 4 pi 



{3/Mi) 



(3/5X? 



(1) 



(2) 



(3) 



where Mi = vfs / \J (5Pi /3pi ) is the Mach number of the 
FS. In the reverse shocked region, by using the strong limit 
jump condition (Landau & Lifshitz 1959; Blandford & Mc- 
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Kee 1976) we have 

^2 _ (r43 + l)[73(r43 



73(2 - 73)(r43 - l)+2 



P3 



73 r4 



73 



P4, 



Ps = (73 - l)(r43 - 1)P3C , 



(4) 



(5) 



(6) 



where r43 — T3T4(1 — P^P^) ~ r4. By using the last remain- 
ing equation ol pressure balance P2 = P3 with 73 = 4/3, we 
can obtain p4 as 



P4 



9 5Mj - 1 



(7) 



20 Mi (r43-i)(4r43+3) 

In order to discuss the physical condition of jets, it is con- 
venient to introduce the dimensionless parameter 



P4 9 5Ml - 1 



Pi 



20 Mi (r43 - i)(4r43 + 3)' 



(8) 



Conventionally we call a jet light for / < 1, and heavy for 
/>!• 

2.1.3 Sideway Expansion 

Most of the hot spot models suppose that the heated plasma 
expands sideway (laterally) to supply the matter in the co- 
coon. Hence, let us discuss the role of sideway expansion 
(escape). 

First, we consider what would occur if there is no lat- 
eral escape. The shock jump conditions shown above do not 
include any information on the longitudinal size of shocked 
region ^3 (see Fig.0. Although little attention has been paid 
on ^3, this strongly reflects a difference between the cases of 
escape and no escape. First, let us show a simple estimate of 
PsRS as a preparation for the estimate of ^3. A jet with rela- 
tivistic speed of r4 decelerates to sufficiently non-relativistic 
speed at the reverse shock, then we have P4RS ~ 1- Besides, 
it is well known that, for a strong shock in a relativistic gas 
in which 73 = 4/3, one has the relation of /33rs x /34RS ~ 1/3 
(e.g.. Kirk & Duffy 1999). Therefore we have Psrs ~ 1/3. 
Next, by exactly solving the Pans by using the relation of 
Pans = (/334 -/3rs4)/(1-/3rs4/334), 0.01 < ft < 0.1, and Eq. 

we obtain the value of Pans is about 0.2 — 0.3. From this, 
we see that accuracy of the simple estimation of P^-rs ~ 1/3 
is good enough. As for a timescale, we take the age of Cygnus 
A as ta,ge ^ lO' yr (Carilli et al. 1998). Then, the scale length 
of reverse shocked region ^3 = tagcW3RS is predicted to be or- 
der of Mpc, which definitely contradicts with the observed 
size. The reason for this is simply because we do not take 
the sideway escape effect into account. 

Next, let us consider the case of including the escape 
effect. In reality, the shocked region expands laterally, too, 
which reduces ^3 significantly. The actual difference will ap- 
pear in the size of the shocked region. In order to get more 
realistic value of ^3, we may use tcsc which expresses the ef- 
fective escaping timescale of shocked plasma from the hot 
spot, because shock accelerated particles do not stay in the 
hot spot during the whole lifetime but escape sideway. Then, 
we must evaluate as 



where tesc is the effective escape time and we believe that 
this is identified as expansion time scale of shocked matter. 
This time scale is estimated by tesc ~ ■^'^ where Rus is 
the hot spot size perpendicular to the jet axis and Vcsc is 
the mean escape velocity of shocked matter. If we take I3 as 
observed size of the hot spot -Rhs, this leads to Vesc of a few 
ten percent of the light speed and shocked plasma remains in 
the spot for only tosc ^ 1/lOOOtagc ^ lO'^yr. This picture is 
matched with the basic concept of synchrotron aging model 
(e.g., Carilli et al 1991) for hot spots. Our next concern 
is the escape velocity. Here we treat the time scale much 
longer than tesc, for the system to be steady. The steady 
state requires that escaping amount of mass and energy from 
a shocked region is the same as those newly injected plasma 
in the region. In this case, the important point is that the 
jump conditions remain the same as pure ID case. Conser- 
vation equations of mass and energy in the sideway direction 
correspond to determine the scale length (surface area) of 
the shocked region and the greatly simplified injection and 
escape balance relations are 



-RhsP3U3RS 

J?HsC^3f3RS 



-Rhs'sPs^csc, 

-RHs'st^S^osc 



(10) 



where U3 is the internal energy density in region 3. An obvi- 
ous but not a unique estimate is I3 -Rhs and Wcsc ~ vans- 
To sum up, on the time scale of order tago, the system can 
be regarded as steady and we expect that it evolves so as 
to keep the pressure and velocity constant and a series of 
the shock jump conditions in the present work remain the 
same ones. Hence we conclude that our ID treatment well 
describes the shock dynamics of FR II sources because the 
shocked plasma escapes in a short time scale. 

On the timescale close to tesc, the study of 2D numerical 
simulation is required to assess the deviation of our simple 
analysis from more detailed hydrodynamical perspective. 

2.2 Electron Acceleration Efficiency 

Here we introduce two important quantities: the ratio of 
mass density of accelerated electrons to that of shocked jet 
(i.e. total particle mass density in the shocked jet) is defined 
as and given by 

P3,acc 



c 



pa 



(11) 



and the ratio of internal energy density of accelerated elec- 
trons to total one of shocked jet is defined as e and given 
by 



<73,acc _ (< 7 > -l)P3,a 
U3 ^ (r43 - 1)P3 



(12) 



^3 ~ tcsct'3RS 



(9) 



where p3,acc and < 7 > are the mass density and aver- 
age Lorentz factor of shock-accelerated electrons, respec- 
tively. The determination of the shock-accelerated electrons 
is shown in the next section. 



2.3 Physical Constraints 

Let us discuss how to give constraints on the energetics 
and plasma composition. To begin with, we estimate the 
minimum Lorentz factor 7min from which value an injec- 
tion process occurs. The determination of 7min is one of 
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the most important issues and it has been discussed from 
various points of view (e.g., Reynolds, Fabian, Celotti & 
Rees 1996; Hirotani et al. 1999; Sikora & Madejski 2000; 
Kino, Takahara & Kusunose 2002 (hereafter KTK02); Asano 
& Takahara 2003). In the case of plasma, the shock 
first converts bulk population of particles into thermal ones 
and accelerates some of them from this thermal pool where 
r^sPeC^ ~ P = UckTc ~ 7minpoC^, then it leads to 7niin — 
r43. On the other hand, in the case of ep plasma, there is 
wide variety of possibilities about 7miii. One extreme case is 
that protons and electrons are separately thermalized where 

r43ppC^ P ~ UpkTp, r43PeC^ ~ Pe = U^kTc ~ 7minPcC^, 

and Te = {mc/mp)Tp, which also leads to 7min — r43. 
The other extreme case is that thermal energy of protons 
is quickly transferred to electrons and they attain one tem- 



perature state r43PpC 



P = 2npkTp 



7minPoC , which 



leads to 7miii — ;;^r43. To make the argument clear, here 
we only treat these three extreme cases as follows; 

(i) The plasma consists of pure (Case I). 

(ii) The plasma consists of pure ep with no energy transfer 
from protons to electrons, which leads to two temperature 
Ti ^ Te (Case II). 

(iii) The plasma consists of pure ep with quick energy 
transfer from protons to electrons and one temperature state 
Ti = Te is realized (Caselll). 

We omit intermediate cases. Summing up, we have 7min as 



7min 



r43 (Case I, Case II) 

^r43 (Case III). 



(13) 



Correspondingly, the rest mass density of shock-accelerated 
electrons is given by 



P3,acc 



??T'C^3,acc 
r'7ma: 



(Case I, Case II) 
(Case III) 



(14) 



where K is the normalization factor of the relativistic elec- 
tron number density and s = 2 is adopted. Using the injec- 
tion rate of relativistic electrons and adiabatic loss time- 
scale icsc, it is written as K — tcscQc (e.g., Mastichiadis & 
Kirk 1997; KTK02). The energy density of shock-accelerated 
electrons is given by 



U3. 



(< 7 > -l)n3,accmeC 



KrUoC In (Case I, Case II) 

(Case III) 



mp 



(15) 



where s = 2 is adopted. 

From the observed non-thermal spectra, we can ob- 
tain the number and energy densities of accelerated elec- 
trons. From the hydrodynamical shock condition, we inde- 
pendently obtain the number and energy densities of total 
particles. Hence we can express e and as a function of a 
single parameter r4. Then the obvious condition 1 ^ is 



rewritten as 
1 ^ 20 



Ml (r43 - 1) n3. 



9 5Mi 



1 Pi 



FS 



3(5Ml-l) Pi V" r43y 
(Case I, Case II) 

4 KmcC^ / 1 _]_\ me 

3(5A1^-1) Pi \ ^43) nip 

(Case III) 

and the condition 1 ^ e can be written as 
20 Ml (< 7 > -1) Jia.acc mc 



(16) 



9 5Ml - 1 

4 KmeC 

3(5M'j-l) 



/3pS 



niCM mp 



In I 



(Case I, Case II) 

mj 



(17) 



3{5M'f-l) Pi 

(Case III) 



■In 



From this, we stress an important point that 1^ and e have 
very weak dependence on r4 through r43. This is understood 
as follows. 

Given the ICM physical quantities and hot spot advance 
speed, P3 is fixed from Eqs. Q and When the jet speed 
is relativistic, we obtain the relation P3 oc r43p3 ~ r4P3. 
Therefore, we have the relation ps oc r4^. On the other 



hand, ps.acc is given by p3,i 



« 7„ 



for s = 2. 



Thus has only a weak dependence on r4. Next we consider 
e. The average Lorentz factor of relativistic electrons is < 
7 >cx 7min ln(7max/7min). Then Combining these relations, 
we can see e cx< 7 > C/r43 oc ln(7niax/r4), therefore e has 
also a weak dependence on r4. 

So far, we have two constraints on the energetics and 
plasma composition of jets, i.e., e ^ 1 and ^ ^ 1. Another 
constraint is 713, tot ^ W3,acc which is written as 



1 > 



W3,i 



n3,tot 



C (Case I, Case III) 

(Case II). 



(18) 



By these three constraints, we investigate the physical con- 
dition of FR II jets in the next section. 



3 APPLICATION 

We explore the energetics related with non-thermal elec- 
trons. As specific FR II sources, we deal with Cygnus A 
and 3C 123. The nearby radio galaxy Cygnus A (z=Q.0562) 
is well known as the very powerful radio galaxy. Recent ob- 
servation with the Chandra X-ray Observatory resolves the 
X-ray emission from the two brighter hot spots (A and D) by 
Wilson et al (2000) with sufficiently high angular resolution 
~ 0.5 . We adopt the observational data of spot A compiled 
by Wilson et al (2000). 3C 123 {z = 0.2177) is famous for its 
high luminosity and peculiar radio structure. The lobes take 
the form of diffuse twisted plumes unlike those in any other 
well-known sources (e.g., Riley & Pooley 1978; Hardcastle et 
al. 1997). We adopt the observational data including Chan- 
dra's one compiled by Hardcastle et al (2001). After this 
preparation, in §3.2, we apply the method in the previous 
section to the representative FR II jet Cygnus A and con- 
strain the energetics and plasma contents. Since the advance 
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speed of hot spot in 3C 123 is not constrained so far, we do 
not treat the case of 3C 123 in this work. We use Ho ~ 50 
km s~^ Mpc~^ and go = throughout this paper. 



3.1 Non-thermal Emission from Hot Spots 

Using the observed spectrum from hot spots, we can deter- 
mine the physical quantities related with non-thermal elec- 
trons. The purpose of this section is to explore the energetics 
in the hot spots based on the multi-frequency spectra. Here 
we identify the hot spot as the region 3 and identify the ICM 
as the region 1, as shown in Figure 

The first point to be discussed is the X-ray emission 
mechanism from hot spots. A thermal model tends to re- 
quire much higher densities than the upper limit obtained by 
the absence of Faraday depolarization in hot spots. There- 
fore, non-thermal emission is most plausible for X-ray emis- 
sion from these hot spots. Next, we must check the rel- 
ative importance of the seed photons for inverse Comp- 
ton scattering. The energy density of locally produced syn- 
chrotron photons in the hot spots is )7syn = f^j^a" = 

8 X 10-" ( io^/;^\-0 (r^)"'«>^g ^ i« the 

size of the hot spot and we set Lsyn = I/syn,o because 
the advance speed of the hot spots is inferred to be non- 
relativistic. The cosmic microwave background (CMB) en- 
ergy density is given by C/cmb = arcMB(l + z)^ = 4.2 x 
10"^'^ (1 -I- z)'*erg cm""^ where a is the radiation constant. 
Compared to L^syn, Ucmb is safely negligible. Alternatively, 
the importance of other seed photon sources for IC scat- 
tering has been suggested for X-rays from radio lobes (not 
a hot spot) (e.g., Brunetti et al. 2001). It is likely that 
there is a strongly beamed emission from the blazar jet 
in the nucleus, which comes along the jet axis and is in- 
jected into the hot spots. Hence it must be checked if this 
process is effective or not. Although we cannot directly ob- 
serve the beamed emission from the inner sub-parsec jet (i.e., 
blazar) of FR 11 radio galaxies, unified schemes for AGNs 
(e.g., Urry & Padovani, 1995) suggest that the inner jets 
should have the same properties as quasar-hosted blazars 
(QHB). The observed bolometric synchrotron luminosity of 
QHB Lsyn, o, core is evaluated by Kubo et al. (1998) and typ- 
ically L: 



syn,o,corc iC^ erg s ^. Thus, the injected syn- 
chrotron photon energy density is )7syn,corc ~ ^=yn,o,<:oro _ 



1 X 



erg cm 



where Dhs is 



, 100 kpc 

the distance of the hot spot from the core. Only when an 
extremely luminous blazar component is hidden and/or Dhs 
is fairly small, (7syn,core component is dominant as seed pho- 
tons. In this paper, we treat the case of L^syn > C4yn,coro 
which seems to be plausible for hot spot sources observed 
by Chandra. (The case of ?7syn < f/gyn.core is discussed 
by Brunetti et al. 2001.) Therefore it seems reasonable to 
suppose that seed photons of inverse Compton X-rays are 
synchrotron-dominated and we adopt the synchrotron-self 
Compton (SSC) model in this paper as in most of the pre- 
vious works. 

Before we get on the issue of SSC model, we describe 
our naming on non-thermal electron distributions used in 
this paper. When the radiative cooling time is shorter than 
the adiabatic expansion loss time of electrons tcsc, the elec- 
tron spectrum steepens by 1 in the power law index because 



the cooling time is inversely proportional to the energy. Ac- 
cording to the value of 7br which describes the break of the 
electron Lorentz factor by radiative cooling, we can classify 
relativistic electron spectra into three regimes. Schematic 
pictures of the electron spectra are shown in Figure |5| If 
7br > 7max, we Call it Weak cooling regime and if jhi < 7max, 
we call it moderate cooling regime, and if 7br < 7min, we call 
it strong cooling regime, in this paper. In comparison with 
the conventional terminology of these cooling regimes which 
Sari, Piran & Narayan (1998) originally named in the study 
of GRB afterglows, strong cooling corresponds to fast cool- 
ing, while moderate cooling corresponds to slow cooling. In 
their study, there is no counterpart for weak cooling. But we 
need to consider this situation for hot spots. Hence, in this 
work, we renamed these cooling regimes. 



3.1.1 Analytic Estimations 

It is useful to estimate the model parameter analytically 
since it gives us a good insight into the physics behind the 
relationship between the model parameters and typical ob- 
servables. Hence we describe some analytic estimates of the 
model parameters following the method of KTK02 and Kino 
(2002) before presenting numerical results. 

In the case of hot spots, we can directly observe the 
emission region size R, and we can set the Doppler fac- 
tor S — 1. The other three parameters, B, 7max, and Qe 
remain to be determined. The three typical observables 
in the observer frame are: v, 



syn,o,max 



the maximum syn- 
chrotron frequency, Lsyn.o total synchrotron luminosity, and 
Lssc.o total SSC luminosity. For example observables of hot 



spots of Cygnus A are Fsj 



1 X 10""ergs"^cm"^, 



and 



Fbbc ~ 5 X 10~'^^erg s~^ cm~^ (obtained from the direct in- 
tegration of the multi-band spectrum). As for the z^syn,o,max, 
we cannot define a strict value from the observed data since 
IR and optical data only give upper limits. Here we tenta- 
tively set fsyn.o.max ~ 10^^ Hz, although some suggest higher 
"cut-off" frequencies close to ~ 10^'' Hz (e.g., Carilli et al. 
1999). Since Eqs. 1141 and 1151 have weak dependence on 
7max, no significant effect is expected on the electron accel- 
eration efficiency. Observed total flux, typical frequencies, 
and luminosity distance are scaled as 

-fsyn^o LgsCjO { 1 f\\ 

10 ^^erg s^ cm lO^'^erg s ^ cm^ 



and 



lOiiHz 



R 
Ikpc ' 



100 Mpc 



(20) 



Using Eqs. (7), (10), and (11) in KTK02 which connect the 
observables and model parameters, we obtain 



B = 5.0 X W /sy„/s 

and 



-1/2 -1 Q 



In 



= 1.3 X 10^ f-'J'f.i>. 



1/2 j-l/4, 1/2 1/2 



(21) 



(22) 



Next, we must discuss the break Lorentz factor 7br which 
appears in the electron spectrum. The one-zone SSC model 
with electron escape on a finite time scale predicts that the 
break Lorentz factor 7br appears at an energy where tcsc = 



tcooi. From this, we obtain 



4((7B+C/3yn)CTT7br 



where Ci 



6 



expresses the adiabatic expansion loss. Then, coohng break 
in the electron spectrum can be written as follows; 



7b. = 1.0 X 10'' /-?Jssc ( 1 + ^ 

/syn 



rci . 



(23) 



where ci — 3Ci and this is reasonable when the escape speed 
is comparable to the relativistic sound speed. In the moder- 
ate cooling regime, we have 



K 



47r_RV3 

„2 



+ ^11 



fssc 
/syn 



-3 + 2 J 
7 07 



-s + 1 



(24) 



where K — qctcsc- For the weak cooling regime, no break 
feature appears in the spectrum. In this case, Qc is written 



K = 



4ttR3/3 



4 

-cttcUb 



-3 + 2 



(25) 

Clearly, if 7br ~ 7max, Eqs. 12411 and 1251 give the same value 
of Qc- 

In order to check which regime is realized for each hot 
spot, we analytically estimate 7br and 7max and compare 
them for ci = 1. Adopted observables are shown in Table 
Qand the resultant 7br and 7max are shown in Table |21 As 
it turns out, weak cooling is plausible for these hot spots 
within this order of estimation. 

By analytic estimation of Eqs. 12311 and 122L we can see 
that weak cooling or marginally moderate cooling is plausi- 
ble for typical hot spots. The case of 7br ~ 10^ - 10^ is not 
appropriate, since the theoretically predicted spectrum for 
this case underlies below the observed radio spectrum. 



3.1.2 Numerical Fitting Results for Cygnus A and 3C123 

Numerical calculations are the best way to obtain the ex- 
act electron energy spectrum and here we show the numeri- 
cal fitting results of observed multi-frequency spectra in hot 
spots. The detailed treatment of synchrotron and inverse 
Compton scattering processes is shown in KTK02. Note 
that, in our model, the injected power law electron spec- 
trum does not have a sharp cut-off but have an exponential 
cut-off such as 7"" exp(— 7/7max) (e.g., Mastichiadis & Kirk 
1997; KTK02). 

The thick solid curve in Figure |3 shows the best fit 
spectrum for the hot spot A in Cygnus A obtained by the 
one-zone SSC model for 7min = 1 and ci = 1. Correspond- 
ing physical parameters are shown in Table |21 Comparing 
these parameter values with those obtained by the analytic 
estimation (Tablej^J, we find that analytic estimation in the 
previous subsection is a good approximation. The important 
result is that f/3,acc/C^B = 7 in the hot spot A for 7min = 1 
(this roughly corresponds to Case I and Case II since is 
an order of unity) . This means that the kinetic power is dom- 
inated in the hot spots similar to the case of TeV blazars 
(KTK02). In order to check the Case III, we examine the 
case of 7min ~ 2000, and this is drawn by a dotted line. In 
this case, we have f/s^acc/C^B = 0.8. As we can see in Figure 



13 predicted fiux is below the observed one below ~ 10^° Hz 
and for Chandra X-ray data. Hence, we rule out the Case 
HI. 

The thick solid curve in Figure |1] shows the best fit 
spectrum for the eastern hot spot in 3C 123 obtained by the 
one-zone SSC model for 7min ~ 1 and ci = 1. Although the 
radio spectrum is rather poorly fitted, this does not much af- 
fect our conclusions as discussed later in this subsection. As 
is for Cygnus A, by comparing Table|31with Tabled we find 
that analytic estimation in the previous subsection is a good 
approximation. The important result is that f/a.acc/C/B = 18 
for 3C 123. This means that the kinetic power is dominated 
in the hot spots similar to the case of TeV blazars (KTK02) 
and Cygnus A hot spot. Similar to the case of Cygnus A, 
we examine the case of 7min ~ 2000, we have Ua^acc/Us ~ 3 
and also we rule out the Case HI for 3C 123. 

To sum up, we find that the energy density of relativistic 
electrons is about one order of magnitude larger than that 
of magnetic field which means that conventional minimum 
energy hypothesis is not strictly true in these hot spots. We 
can also find that relatively low 7min ~ a few is preferred 
based on the SSC model fitting. In this case, energy den- 
sity of non-thermal electrons is about one order of magni- 
tude larger than that of magnetic fields. Hence, we conclude 
that the assumption of neglecting the magnetic effects in 
the shock jump condition in §2 is a correct one. Additional 
new finding in this spectral calculation is that the twice- 
scattered Compton component is predicted for these sources 
in the GeV energy band. If this bump can be detected in 
future observations, it must give strong constraints on phys- 
ical quantities in these hot spots. However the predicted fiux 
does not reach the detection threshold of GLAST. 

Compared with previous works, for Cygnus A, the re- 
sultant physical quantities and energy densities in this paper 
are almost the same as the result of Wilson, Young, & Shop- 
bell 2000) . According to their adopted value for the hot spot 
A, their result leads to (73,acc/f/B ~ 4. A difference between 
the present work and the work of Wilson, Young, & Shopbell 
2000 is that they adopted 7br — 4 x 10'' , while we show that 
weak cooling regime is consistent with observations (see be- 
low as for a solution in moderate cooling regime) by solving 
the electron kinetic equation both analytically and numer- 
ically. We find that low radiative efficiency by comparing 
relativistic electron energy and radiation field energy. For 
3C 123, Hardcastle, Birkinshaw, & Worrall (2001) regard 
that magnetic fields strength in the hot spot is consistent 
with being equipartition by assuming 7min — 1000. This 
roughly corresponds to the case of ep plasma composition 
with one temperature for which 7min ~ mp/mo. This as- 
sumption seems to contradict with their another assumption 
that no relativistic protons are contained in the spot. 

Let us examine whether the above results depend on 
the assumption of ci = 1 since predicted radio spectra show 
only a global matching with the observed one. Some of the 
previous work reported the existence of the spectral break in 
the radio band for Cygnus A (Meisenheimer et al. 1997) and 
3C123 (Looney & Hardcastle 2000). To reproduce a cooling 
break (to enter the moderate cooling regime), the smaller 
value of ci is required. Hence we examine the case of smaller 
ci. We select the case of ci = 1/6 which corresponds to 
the case of escape velocity 0.055c for the spot A (Carilli 
et al. 1999). Although most of the effects of smaller ci are 
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absorbed by a decrease of go, the decrease of high energy 
electrons leads to a decrease of high energy portion of the 
synchrotron spectrum. To compensate for this, somewhat 
higher 7max and B will be required compared to the case of 
ci = 1. The numerical results are tabulated in Tables 3 and 
4 and predicted spectra are shown in Figs. 3 and 4 by thin 
solid curves. 

For Cygnus A, the differences in the resultant parame- 
ters between the cases of ci = 1 and ci — 1/6 are very small 
as is seen in Table 3. The spectral fitting reproduces equally 
well the observations for ci = 1/6 and 7min = 1 as shown in 
Fig. 3 by the thin solid line. The case with large 7min under- 
predicts low frequency radio and SSC X-ray fluxes as is seen 
in Fig. 3 by the thin dotted line and is ruled out. For 3C123, 
the differences are within a factor of 2 as is seen in Table 4; 
the magnetic field strength is 45% larger and f/a.acc/f^B is 
smaller by a factor of 2 for ci = 1/6 than for ci = 1. How- 
ever, the spectral fitting is not significantly improved even 
for ci — 1/6 and 7min = 1 as shown in Fig. 4 by the thin 
solid line. Although the prediction better reproduces high 
energy radio spectrum, it is well above low frequency radio 
data. For ci = 1/6 and 7min = 2000 whole synchrotron spec- 
tra are better fitted but the SSC fiux in the X-ray band is 
below the observation. In this sense, ci = 1/6 and interme- 
diate 7min may best reproduce the observed data, although 
we do not further pursue this point in this paper. If this is 
the case, U^^^cc/Ub becomes smaller, but still greater than 
1. 



3.2 Physical Condition of Cygnus A 

3.2.1 Known Quantities 

Here we investigate the plasma content and electron ac- 
celeration efficiency in Cygnus A hot spot which is well 
known as the best candidate for studying the nature of FR 
II radio source. Regarding the unshocked ICM, the con- 
tinuous study of the X-ray observations reveals its tem- 
perature and number density (Arnaud et al. 1984; Ueno 
et al. 1994; Smith et al. 2002). Here, we adopted the the 
number density and temperature of Cygnus A cluster as 
pi = nipnicM — 10~^mp g cm~^ and Ti ~ 4 x lO'' K shown 
in Arnaud et al. (1984). Then we have the ICM pressure by 
Pi = 2nicMfcricM. The velocity of the fluid in the upstream 
is 111 = by definition. Hence we can obtain three down- 
stream quantities of the forward shock based on the ICM 
observations. About the unshocked jet, we assume P4 = 0. 
The advance speed of the hot spot /3hs has been estimated 
by synchrotron spectral aging methods (Carilli & Barthel. 
1996; for review), 02 = Pd. = /3hs- We use the value of «3 
taking into account an uncertainty of about one order of 
magnitude. In the present work, we examine the case of 
0.01c < U3 < 0.1c. In other words, we take account the 
uncertainty in 113 instead of B as is the synchrotron aging 
methods (e.g., Carilli et al. 1991). The velocity of the un- 
shocked jet TiPi is not directly measured. The upper limit of 
the jet velocity is inferred from the initially ejected jet speed 
corresponding to the sub-parsec scale jet (i.e., blazars) ve- 
locity and the upper limit of r4 is taken to be about 30 
(e.g., KTK02). With regard to the lower limit of the jet 
velocity, we set the mimimum value as r4 > 3 for simplic- 
ity. Note that recent statistical studies of measurements of 



jet-to-counter jet flux ratio in FR II sources imply that the 
minimum r4 is close to unity (e.g., Wardle & Aaron 1997; 
Hardcastle et al. 1999). Thus we can fix 5 quantities such as 
pi. Pi, ui, vz and P4. Then, we solve 6 physical quantities 
such as p2, P2 = P3, P3, P4, wfs, and wrs as a function of 
r4 by following the procedure shown in the previous sec- 
tion. As for the non-thermal electrons, based on the result 
in the previous subsection, we examine the case of s = 2.1, 
7max = 1 X 10*, and K = I x 10~^ cm"^ based on Table 
3. As for the relativistic electron number density 713, acc and 
normalization K, here we use the value of 7min = 1 in Table 
3, then we have ns^acc ~ ^7min- O^'^ resultant value of K 
and the one obtained by Wilson et al. (2000) is consistent 
with each other. 



3.2.2 Mass Density of the Jet 

Figure |5] shows rest mass density ratio of the unshocked 
jet to the unshocked ICM. It should be stressed that p4 is 
not a free parameter but a solvable quantity as is shown 
in §2. For r4 ^ 1, one can see p4 oc Pyis^az- ^iS- IHl 
the shaded region corresponds to 0.01 < /3hs < 0.1. At 
r4 ~ 3, /3hs = 0.1 (which corresponds to relatively large 
Mach number Mi ~ 50) leads to / ~ 10"^ while /3hs = 
0.01 (corresponds to Mi — 5) shows / ~ lO^''. Given r4, 
allowable range of the rest mass density p4 spans about two 
order of magnitude because /3hs is uncertain by one order 
of magnitude. For fixed /3hs, as r4 increases, the density 
decreases according to p4 cx V^^ oc V'^^ . From this, we can 
directly see that the jet is very light for all range of r4. 

Let us compare these results with the previous works. 
Many authors have done the numerical simulation of extra- 
galactic jets (e.g., Norman, Smarr, Winkler & Smith 1982; 
Clarke, Norman, & Burns 1986; Marti et al. 1997) and it 
has been found that light jets are required for producing 
cocoon structures (e.g., Norman, Smarr, Winkler & Smith 
1982; Komissarov & Falle 1998). Although most of previ- 
ous works have studied mainly the range of 10~^ < / (e.g., 
Cioffi & Blondin 1992), a study of very light jet has been 
recently reported by Krause (2003). He examined the range 
of lO"'^ < / < 10"^ which has not been studied so far but 
matches with our result. Based on the morphology of Cygnus 
A obtained by Chandra (Smith et al. 2002), they claimed 
that / < 10^^ which agrees well with our result although 
their work and ours are definitely different approaches for 
estimating the mass density of the jet. We should add to 
note that, the lightness of the jet is consistent with the rel- 
ativistic speed of the jet and the non-relativistic advance 
speed of the hot spot. 

3.2.3 Total Kinetic Power of the Jet 

Figure |S| shows the jet power compared with the Eddington 
power. The total jet kinetic power is defined as 

Lkin = 7rPHScr|/34P4C^ (26) 

Similar to Fig. |H1 the important point is that total kinetic 
power Lkin in Cygnus A is not a free parameter but consis- 
tently obtained by solving p4 for given r4. As Lkin oc 0^^, 
the allowed range of the total jet kinetic power spans about 
two order of magnitude. The reason for weak dependence 
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on r4 is already shown in §2, so we do not repeat it here. 
Compared with the Eddington power, when a hot spot has 
relatively high advance speed of 0.1c, the jet power exceeds 
the Eddington power even for a relatively heavy black hole 
mass 10® Mq. This may suggest that the hot spot advance 
speed is slower than 0.1c. 

Compared with the other FR II sources, kinetic power 
of Cygnus A with Phs = 0.1 is by a factor of a few larger 
than that of brightest FR II sources reported in Figure 1 of 
Rawlings & Saunders (1991). So Cygnus A is one of the most 
brightest FR II sources in low-redshifts z < 0.5 not only 
in the radio band (Carrili & Barthel 1996) but also in the 
kinetic power. Related to this, flat spectrum radio quasars 
(FSRQs) are believed to be the inner core part jet (typically 
sub-pc scales) of the FR II sources (Urry & Padovani 1995) . 
Then, the total kinetic power of this core region estimated 
by the SSC analysis as Lkin^corc ~ isyn,corc,o/r^ ~ 10*^ erg 
s~^. Here we used the data of brightest FSRQs Lsyn.corc.o ~ 
10*^ erg s~^ in Kubo et al. (1998) and assume F ~ 10. Using 
this, 

10-' < ^±pS2Il < 10-1 (27) 

Within the framework of internal-external shock scenario 
(e.g., Piran 1999), the value almost equals to the 

ratio of the dissipation rate by the internal shock to that 
of the external shock, only a small fraction of total kinetic 
power from the central engine is dissipated in sub-pc scales 
and most of the bulk kinetic power remains up to lOOkpc 
scales and then is dissipated by strong deceleration by the 
ICM pressure. If the adopted core luminosity estimate is 
acceptable, then our estimate shows suitable agreement with 
the internal-external shock scenario similar to the case of 
GRBs. 

3.2.4 Electron Acceleration Efficiency 

Figure |7| shows the energy density ratio of accelerated elec- 
trons to total particles in Cygnus A hot spot (spot A). In 
the same way as mass density and kinetic power of the jet, 
allowable range of e spreads over two order of magnitude as 
about 10"^ < e < 1 because of the uncertainty of hot spot 
advance speed. We see that the case of /3hs < 0.01 does not 
satisfy the constraint Eq. (1171 and is ruled out. It should 
be noted that if we definitely determine /3hs of Cygnus A 
as close as 0.01 in the future, then we must reconsider the 
shock structure taking the back-reaction of accelerated elec- 
trons into account (e.g., Drury & Voelk 1981; Malkov 1997). 
If /3hs 0.1 is confirmed, then we can employ the simplest 
test particle theory for the Cygnus A spot. It is also worth 
noting that in the case of the rapidly advancing hot spot 
Pi ~ 0.1 the result shows the existence of missing thermal 
power which is about two orders of magnitude larger than 
that of accelerated electron kinetic power. 

It is well known that there is a major gap which sep- 
arates numerical simulations from direct comparison with 
observed image. The gap is that we have little knowledge 
about the relativistic electron production rate (e.g., Krolik 
1999). We again stress that one of most important achieve- 
ments of our work is that we show the simple method to 
remove this difficulty and estimate electron acceleration ef- 
ficiency. 



3.2.5 Plasma Composition 

Figure |H| shows the number density ratio of accelerated elec- 
trons to total particles in the hot spot. In the case of 
plasma, we estimate ns^acc/ji-s = C- In the case of two- 
temperature ep plasma, we estimate ns^acc/wa = 1837(^. In 
the spot in Cygnus A, we find that if the plasma is com- 
posed of a two temperature ep plasma (Case II), the num- 
ber density of accelerated electrons exceeds that of ther- 
mal particles. Fermi acceleration predicts that some frac- 
tion of thermal particles is converted into relativistic par- 
ticles. Hence we can exclude the case of ns^acc/wa > 1 by 
definition. Hence, we can derive the important result that 
a two temperature ep plasma composition is ruled out. In 
Figure 1^ we summarize this plasma composition diagnosis. 
Combined with the result of rejection of one-temperature ep 
plasma (Case III) by the spectrum fitting result for Cygnus 
A, we draw the conclusion that only plasma is acceptable 
for the Cygnus A hot spot. 



4 CONCLUSIONS 

To sum up, we investigate the energetics and plasma com- 
position in FR II sources using a new simple method of 
combining shock dynamics and radiation spectrum. For sim- 
plicity, we examine three extreme cases of pure electron- 
positron pair plasma (Case I), pure electron-proton plasma 
with separate thermalization (Case II), and pure electron- 
proton plasma in thermal-equilibrium (Case HI). 

Based on the SSC model with escape velocity as ci = 
3ucsc/c = 1 (roughly corresponds to weak cooling) and ci = 
1/6 (roughly corresponds to moderate cooling), we estimate 
the number and energy densities of non-thermal electrons 
using the multi-frequency radiation spectrum of hot spots. 
The results of SSC analysis for the hot spots of Cygnus A 
and 3C 123 (see Figures 01 and m Tables 2 and 3) are as 
follows: 

• The energy density of relativistic electrons is about 10 
times larger than that of magnetic fields. 

• For Cygnus A, we find that both of the ci = 1 and 
ci = 1/6 with 7min = 1 explain the observation. In both 
cases, the radiative efficiencies are low (typically ~ 1%) and 
resultant parameters are very similar. Case HI is not accept- 
able because predicted photon spectra do not give a good 
fit to the observed one. 

• For 3C123, ci = 1 with 7min = 1 gives a global fit to 
the observations although detailed radio spectrum is rather 
poorly fit. The ci = 1/6 of 7min = 2000 gives a better fit 
on the radio spectrum but underpredict the X-ray fiux. The 
difference between Cygnus A and 3C123 is an interesting 
open issue. 

Independently, with the ID shock jump conditions taking 
account of the finite pressure of hot ICM, we estimate the 
rest mass and energy densities of the sum of thermal and 
non-thermal particles in hot spots. We utilize the condition 
that the obtained rest mass, internal energy, and number 
densities of non-thermal electrons should be lower than those 
of the total particles determined by shock dynamics. The 
results of the energetics and plasma composition by using 
our new method for Cygnus A in the range of 3 < r4 < 30 
are as follows: 
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• We estimate the ratio of jet rest mass density to that 
of ICM by solving the shock jump conditions. We find that 
the jet is very light 10"'^ < / < 10"^. 

• The total jet kinetic power is examined. We found that 
it has the kinetic power of about 10''® < Lkin < 10''^ erg s~^. 

• We examine the electron acceleration efficiency in the 
hot spot. We found that the efficiency is about 0.01 < e < 1. 

• Plasma composition is investigated. We find that both 
pure two temperature ep plasma (Case II) and pure one 
temperature ep plasma (Case III) can be ruled out. Hence 
we conclude that the plasma composition is most likely to 
be dominated. 
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Table 4. Physical paraiiietorfs from SSC analysis for -30123 



Notes: Fgyn and Fssc are the bolometric flux. 

Table 2. Analytic estimate of physical parameters in the hot 
spots. 
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Table 3. Physical parameters from SSC analysis for Cygnus A 
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< 7 > is the average Lorentz factor of non-thermal electrons and 

7min = 1. 



Parameter 


o01z3(Ci = i) 


oUlzo(Ci = 1/6) 


B{G) 

7m ax 

qe(cm-3s-l) 


1.0 X 10-* 
5.0 X 10* 

2.1 X 10-15 


1.45 X 10-* 
8.0 X 10* 
4.4 X 10-16 


n^3,acc( cm-3) 
< 7 > 


1.3 X 10-3 
6.5 


1.6 X 10-3 
6.1 


-P"syn(crg S-l cm-2) 
Fssc(erg s-l cm-^) 
-^-o,kin{erg s-l) 
-'^poy{erg s-l) 

-'^rad(crg s-l) 


2.0 X 10-12 

2.8 X 10-13 
5.8 X 10*6 
3.2 X 10*5 
5.1 X 10** 


2.4 X 10-12 
1.9 X 10-13 

6.6 X 10*6 
6.8 X 10*5 

5.7 X 10** 


r/rad/f/3.acc 
(:!.„■.■/' B 


8.9 X 10-3 
IS 


8.7 X 10-3 
10 



Notes: -Fgyn and Fssc arc the bolometric flux. 

< 7 > is the average Lorentz factor of non-thermal electrons and 

7inin = 1. 
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Figure 1. Schematic view of shock dissipation at the head region 
of the relativistic jet. Upper panel shows the model geometry of 
shocks in an FR II radio galaxy. Lower panel shows the mass 
densities in each region. The reverse shocked region is identified 
as a hot spot. The forward shock is identified as a bow shock. The 
cocoon is composed of shocked jet material that has expanded 
sideways. The length I3 is the longitudinal size of the hot spot. 
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Figure 2. Cooling regimes for non-thermal electrons. For strong 
cooling (so-called fast cooling in Sari, Piran & Narayan 1998) , in- 
jected electrons immediately lose their energy by radiative cool- 
ing. For moderate cooling (slow cooling), at high Lorentz factors, 
radiative cooling decreases the number density of electrons and 
leads to a break in the spectrum. For weak cooling, the cooling 
time scale is sufficiently long and the radiative cooling docs not 
change the injected electron spectrum. An example of specific 
objects in each regime is also shown. 




Figure 3. One-zone model spectrum for the emission from the 
hot spot A of Cygnus A. The observed data are taken from those 
compiled by Wilson, Young, & Shopbell (2000). The thick solid 



line shows the best fit spectrum in the case of 7^, 



= 1 and 



ci = 1. The thick dotted line shows the case of 7inin = 2000 and 
ci = 1, which is not acceptable to fit the observed one. The thin 
solid line shows the case of 7min = 1 and ci = 1/6. The thin 
dotted line shows the case of 7min = 2000 and ci = 1/6 can be 
also ruled out. 
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Figure 4. One-zone model spectrum for the emission from the 
hot spot in 3C 123. The observed data are taJcen from those com- 
piled by Hardcastle, Birkinshaw, & Worrall (2001). The thick 
solid line shows the best fit spectrum. The thick dotted line shows 

the case of 7min = 2000, which is not acceptable. The thin solid 
line shows the case of 7niin = 1 and ci = 1/6 which overestimate 
the observed radio flux. The thin dotted line shows the case of 
7min = 2000 and ci = 1/6 which exlpains the radio band while it 
tends to underestimate the X-ray band. 



Mass Density of Cygnus A Jet 

(n,= 10"'cm"^T,=4*10' K) 




5 10 15 20 25 

(Lorentz factor of the jet) 



Figure 5. The rest mass density ratio of the AGN jet to the hot 
ICM in the case of Cygnus A (see also Figure 1). Proton number 
density and temperature in the ICM (region 1) are taJcen from 
Arnaud et al. (1984). 
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' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 
Jet Power of Cygnus A 

(n, = 10"^cm"^Tl=4*10'K,R=^ kpc) 

^3=0.1 U,=54) 
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n3..cc/n3.ioi in Cygnus A (Spot A) 
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Figure 6. The total kinetic power of Cygnus A jet. Sliadowed 
region corresponds to tlie allowable range of 0.01 < P3 < 0.1. 
To generate higher hot spot speed /Ss, higher pi is required. We 
show the Eddington luminosity (dotted lines) for comparison for 
several values of black hole mass. 



Figure 8. The number density ratio of ns.acc/ns.tot. The upper 
shadowed region corresponds to the range of 0.01 < /Ss < 0.1 for 
pure ep plasma with two temparature (Case II). This region is 
ruled out simply because the value fi3,acc/^3,tot > 1 is forbidden 
by definition. The lower shadowed region corresponds to the range 
of 0.01 < (83 < 0.1 for purely plasma (Case I) and pure ep 
plasma with one temperature (Case III) . 
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U3..co/U3,t„i in Cygnus A (Spot A) 
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Figure 7. The energy density ratio of accelerated electrons to 
the total particles. Solid lines represent the ratio in the case of 
= 0.1 and (83 = 0.01, respectively. 
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Figure 9. The summary of plasma composition diagnosis and 
sketches of the electron energy spectra in Case I, II, and III. 
Thick solid line in each panel shows the spectra in the hot spot. 
Thin solid line in each panel shows the cold jet. Dashed lines 
in Case II (in left lower panel) imply that the number density 
of non-thermal electrons exceeds that of total particles and that 
this case is ruled out. The normalization of non-thermal electron 
number density in the hot spot is fixed by SSC model fitting. 



